Europhysics Letters 



PREPRINT 



Is the entropy at the liquid-gas critical point of pure fluids 
proportional to a master dimensionless constant? 

Yves Garrabos 

Equipe du Supercritique pour I'Environnement, les Materiaux et I'Espace - Institut de 
Chimie de la Matiere Condensee de Bordeaux - Centre National de la Recherche Scien- 
tifique - Universite Bordeaux I - 87, avenue du Docteur Schweitzer, F 33608 PESSAC 
Cedex France. 



PACS. 64.60.-i - General studies of phase transitions. 
PACS. 05.70. Jk - Critical point phenomena. 
PACS. 64. 70. Fx - Liquid- vapor transitions. 



Abstract. - From a minimal set made of four scale factors defined at the liquid-gas critical 
point of a pure fluid, and one adjustable parameter which accounts for particle quantum effects, 
we demonstrate here a master singular behavior of the correlation length for the one-component 
fluid subclass, using an asymptotic scale dilatation of the physical fields. Such master behavior 
observed within the preasymtotic domain is in conformity with the renormalized $^=3 Field 
Theory predictions at large correlation length scale of the fluctuating order parameter, for the 
complete universality class of the symmetrical uniaxial 3D-Ising-like systems. The following 
consequences are discussed: (i) A comparison between the critical state of pure fluids and the 
zero-temperature state leads to an intuitive analogy with the (Nerst) third law of thermo- 
dynamics, which authorizes specific master form for hyperscaling within the subclass of pure 
fluids; (ii) A master constant value of the non-dimensional critical entropy can exist for all 
the pure fluids at the short-ranged lengthscale of the molecular interaction. From this latter 
hypothesis, we show that the needed four scale factors are the four preferred directions ex- 
pressing complete thermodynamic (linear) continuity crossing the liquid-gas critical point on 
the (pressure, volume, temperature) phase surface. 

The liquid-gas critical point (CP) of a pure fluid is an unattainable single point on the 
p,Vp,T (pressure, particle volume, temperature) phase surface, in which an infinite degrees 
of freedom are coupled. As a matter of fact, the critical state at p c , T c , and n c — ^f- 
(the critical density number) of a critical amount of matter N = N c filling the volume V, 
is a thermodynamically unstable state, due to the diverging character of the spontaneous 
fluctuations of extensive variables (p is dual to V) . This critical state appears then as a limit 
of thermodynamic stability at which all the stability determinants (the second derivatives 
OU dn )■ th &t were strictly positive for any point of the four-dimensional (4-D) characteristic 
surface 4>u (f, f2j) = 0, become zero for CP (see for example [1]). U (Qj) is the total internal 
energy, while f2j = (S, V, N) are the three associated natural variables. S is the total entropy. 

A comprehensive understanding of the diverging character of the spontaneous fluctuations 
of extensive variables, which characterizes the critical behavior of the one-component fluid 
close to its unstable critical state, comes from the field theory (FT) framework (see for example 
[2]). This theoretical scheme accounts for the infinite degrees of freedom throughout the 
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Hamiltonian of the so-called $^_ 3 (n = l)-model of the symmetrical uniaxial 3D-Ising-like 
systems, with associated coupling constant M4 > and finite cutoff wave number Ao. The 
selected pair of relevant scaling fields, made of the thermal field t weakly fluctuating and the 
ordering field h more strongly fluctuating, takes exact zero- value at the isolated non-Gaussian 
(Wilson-Fisher) fixed point {t = 0; h = 0}. At this fixed point, the fluctuations of the (scalar) 

order-parameter (OP) m = ($), conjugated to h, can reach infinite size ~ L = (V)' 3 ~ 00 j. 

This hypothetical situation at T — T c is then only expected for an unphysical pure fluid 
of infinite size (V" — ► 00) and infinite number of particles (N — > 00), but with finite number 
density n = y = (vp)^ 1 , such that n c — ( y ) cp ^ t^7" Asymptotically close to this 
critical point, i. e. for small values of t and h, power law behaviors of universal features 
are expected whatever the selected system. So that, the set of system-dependent parameters 
which characterizes each one-component fluid is made from i) the critical parameters, such as 
T c ,p c , and n c , ii) the inverse cutoff wave number which characterizes a discrete struture of the 
fluid particles with spacing (Ao) 1 , and iii) the two-scale factors which relate analytically t 
and h to the respective physical fields proper to each system [3] . Starting from this asymptotic 
description of the two-scale universality in the close vicinity of CP, we have postulated in [4] , 
that the characteristic set of fluid-dependent parameters corresponds to the minimal set of 
measured critical parameters needed to localize CP on the normalized (i.e. 3D) p, vp, T phase 
surface. Here normalization refer to particle properties for standard thermodynamics written 
for a constant amount N of matter. Our particle notation, i.e. vp = jj, uses small letters 

with explicit subscript p. The minimal set reads Q™ m = } T c ; up^ c ;p c ; 7 C = ( gfu r [4,5]. 

I ^ ' Vp.c J CP 

We can then make dimensionless the thermodynamic and correlation functions of any one- 

1 

component fluid, using the scale factors, (/3 C ) _1 = k B T c for energy unit, and a c = ( kE p Tc ^j 3 
for length unit (ks is the Boltzmann constant). a c is not dependent of the container size 
L = (V) d [6] and takes a clear physical meaning: a c is the spatial extent of the short-ranged 
(Lennard-Jones like) molecular interaction. Therefore, v c ,i = Tf , is the microscopic volume 
of the critical interaction cell (CIC). Introducing the two characteristic dimensionless number, 
Z c = , and Y c = 7 c fr- — 1, leads to rewrite the minimal set in the more convenient form, 
qmin = j^. ^ Zc; Yc } . Now, (Z^ 1 = n c v cJ is the number of particles that fill the CIC. Z c 
and Y c can then take their physical meaning of two-scale factors to formulate the dimensionless 
master behavior of all the one component fluids asymptotically to their CP [5,7]. 
"Master" singular behavior of the one-component fluid subclass 
As a matter of fact, asymptotic master singular behavior of dimensionless potentials and 
dimensionless size of OP fluctuations only occurs when an appropriate scale dilatation method 
(SDM) is applied to the physical fields, At* = k B l3 c (T -T c ), Ah* = (3 C (fj, p - fj, PiC ), and 
Am* = {a c ) d (n — n c ). fip, dual to N, is the chemical potential per particle. In SDM, the 
"renormalized" fields T*^, and are proportional to the physical fields, in complete 

analogy to the FT framework near the Wilson-Fisher fixed point [3]. For the complete one- 
component fluid subclass, these renormalized fields are now [8] defined by 

T* f = T* = Y c At* 

H* qf = (A* e ) 2 7Y* = (A*,) 2 (Z c )-* Ah* (1) 
M* qf = A* qe M* = A* e (Z c )iAm* 

In Eqs. (1), A* e = 1 + A c accounts for quantum effects on the cut-off parameter for 
T = T c [8]. Writing A c = A„ / (—), [A Tc = ^ r is the thermal wavelength at 

V a " I ' (2Trm r k B T c )2 
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T = T c and hp is the Planck constant], leads to a relative quantum correction of the range of 
molecular interaction at CP proportional to the ratio ^p^- Then X q j appears as an adjustable 
numerical prefactor which incorporates the quantum particle statistics. 

When critical properties of one standard fluid (xenon [9]) are known, Eqs. Q permit 
to define the constant amplitude values of the master critical behavior [5] in the so-called 
preasymptotic domain (PAD) [10], where the singular power laws expressed at the first-order 
of the Wegner expansion [11], are expected to be valid. Conversely, SDM is then able to 
estimate all the critical amplitudes appearing in the [two-terms] Wegner expansions for any 
pure fluid [5,9]. 
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Figure 1 - a) Log-Log scale of £ + (in urn) as a function of T — T c > (in K), along the critical 
isochore, for Xe, Kr, Ar, CO2, SFq, D2O, and 3 He (see the respective colors in the inserted Table); 
b) dimensionless behaviors in units of a c and (/3 C ) _1 , showing failure of the classical corresponding 
state scheme; c) master dimensionless behavior of the renormalized correlation length , as a 
function of the dilated thermal field T* [see Eq. (2)]. The arrow indicates the order of magnitude 
of the expected extension of the preasymptotic domain. Each fluid length scale a c is given in the 
inserted Table. 

To demonstrate this important feature where any leading and first confluent amplitudes 
can be estimated only using Q™ m and A* e , we are here concerned by the critical behavior 
of the actual correlation length £ + along the critical isochore above T c [12]. When the fluid 
cutoff wave number Ao is such that AoA* e = ^- [8], the renormalized correlation length writes 

l* q f = Ao£ = -jfcr = jt~. Within PAD, the two-terms master divergence of £*'^~ reads as 



' gf 
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where v = 0.6304 ± 0.0013 and A = 0.502 ± 0.002 are universal critical exponents [13]. 



The leading amplitude = 0.57 and the first confluent amplitude Z t 



0.385 have 



constant values for the pure fluid subclass (see [12] for detailed analysis and the Refs. [5,7,9] 



for amplitude values of standard critical xenon). Here, postulating Z t 



(i),+ 



cte for the pure 



4 



EUROPHYSICS LETTERS 



fluid subclass, leads to consider the simpler situation in the $^_ 3 (n — l)-theory, where the 
starting point for U4 > (in usual renormalized trajectories), is certainly very close to the 
ideal trajectory between the Gaussian and the Wilson-Fisher fixed points [14]. 

The published fitting results of the correlation length measurements of Xe, Kr, Ar, CO2, 
SF 6 , D 2 0, and 3 He, [15] have been reported on Figure la using dimensioned quantities, which 
make clearly distinguishable each specific fluid behavior. Figure lb gives a representation of 
the dimensionless quantities obtained from the classical theory of corresponding states, using 
(/3 C ) 1 and a c units. The failure of the classical theory is evidenced from the importance 
of quantum effects in 3 He and molecular interaction effects in D 2 0, when comparison to 
the standard monoatomic Xe is made. Final representation of the master behavior obtained 
from SDM is given on Figure lc. The scatter between the data corresponds to the estimated 
precision (10%) on the determination of each fluid correlation length. 

Since a c is a measure of the mean range of interaction forces, Eq. (J2J also provides an 
easy control of the effective extension of the critical domain where the mandatory condition 
£*'j~ ~> 1 is expected to be valid. Then Figure ^ indicates also the order of magnitude of 
the PAD length (at least up to T* < 0.01) where Eq. @ is valid (see also [7,12]). Such 
result confirms previous similar conclusions based on careful analyses of crossover models and 
generalized critical fluid e.o.s. [5,9,16]. 

Characterization of the thermodynamic "CP vicinity" 

To reach CP from a thermodynamic approach needs to perform an infinite number of 
transformations between infinite number of near-critical equilibrium states, leading to infinite 
time to obtain the critical state. By analogy with the Nernst principle (the so-called third 
law of thermodynamics), we can also reformulate the previous sentence as follows. It is 
impossible by any procedure, no matter how idealized, to reach exact critical state of any pure 
fluid in a finite number of operations. In a finite number of operations (or finite time), it is 
then only possible to border the critical point as close as possible, and the final near critical 
state at incipient equilibrium is imposed by the finite sized container. For pure fluids in 
absence of external field, this finite size of the container is measured by i) finite extensive 
values of V and N, fixing the (small) mean value of n — n c (or Vp — Vp, c , equivalently) , 
proportional to the physical OP; ii) the finite intensive value of T, that fixes the (small) 
mean value of T — T c , proportional to the independent physical thermal field. Therefore, the 
finite intensive value of [ip — jj,p tC (or p — p c , equivalently), proportional to the independent 
physical ordering field, is also fixed (from thermodynamics principles). This finite critical 
size of the container governs the natural way for the unstable critical fluid to reach incipient 
stability on an equilibrium state acted by this near critical container, the so-called reservoir in 
statistical mechanics. For such a thermodynamic equilibrium very close to the CP, all the total 
characteristic potentials II = {U, H, A, G, J} [17] are homogeneous functions of the first order 
in terms of their three natural variables among Q< = {S, V, N, T,p, fip}, and Euler's thoerem 
applies [1]. Correlatively, the Gibbs-Duhem equation, SdT — Vdp + Nd^ip — 0, requires 
a normalized description leading to the 3-D representation of thermodynamic equilibrium 
states. As N and V are two independent extensive variables, in addition to the standard 
normalization per particle mentionned above, another equivalent normalized scheme occurs 
using densities for a system at constant volume. The density notation, i.e. n — y or n = y, 
uses small letters. 

When the above first-order (equilibrium) scheme is used in the "CP vicinity", the finite 
critical parameters are basically constitued from the non-zero values of appropriate "first 
order" derivatives. On this thermodynamic point of view, all the critical parameters, the 
so-called preferred critical directions in the following, reflect topological continuity of thermo- 
dynamics at the CP, directly associated to the proper analytic continuity of the three intensive 
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variables T, p, and /ip when "crossing CP". For example of our main present concern, the 

normalized value of Sp^ c = [or s c = (ff 1 ) L should be defined as one among the 

preferred directions. Obviously, the extended set made with all the critical directions (see be- 
low) , includes our minimal set hypothetized as containing the needed information to describe 
critical phenomena (i.e. the needed information to calculate all the leading and first confluent 
amplitudes of the power law singularities of "second order" derivatives). This raises the fol- 
lowing question resulting from thermodynamic equivalence between unattainable critical state 
and unattainable absolute zero state: Is the entropy at the exact liquid-gas critical point of 
any pure fluid a characteristic value of each fluid particle, related to a master dimensionless 
constant reflecting universal fluid nature approaching asymptotically CP? 

Close to the critical point, we are mainly concerned by the fluctuations of N (or V, 
alternatively) . When the energy and the number of particles of a fluid in contact with a particle 
reservoir can fluctuate, the basic link between statistical mechanics and thermodynamics of a 
fixed volume of matter, where T and up are the operating variables can be correctly established 
only from the Grand canonical statistical distribution. The Grand potential J (T, V, lip) = 
—p V is then naturally selected to characterize the equilibrium state of the system maintained 
at constant volume. Its equilibrium state is characterized by the Grand potential density 
j (T,fip) = y, whose opposite, i.e. p(T,fip), gives the characteristic surface close to the CP 
schematized in figure 2 [18]. This result complements the Canonical statistical description 
connected to the Helmholtz free energy A(T,V, N), where the Helmholtz free energy per 
particle dp- (T, Vp) = characterizes the equilibrium state of the system of constant amount 
of matter in contact with an energy reservoir (a thermostat) . To discuss differences between 
minimal and extended sets, needs to consider the e.o.s. pairs, {p (T, Vp) ; s p (T, Vp)}, associated 
to dp, or, {n (T, up) ; s (T, up)}, associated to j(T,fj,p). However, only their common Hp-;T 
diagram contains three new unmeasured characteristic parameters, in addition to the measured 
ones in the usual p; T diagram (see the respective binary diagrams constructed in Figure 2): 

i) Mp, c = 9p, c = Pc^p, c , with a,p. c = jp. c + gp, c = -p c Vp, c + Hp.c and a± c = j± <c + g^ c = 
—Z c + fj£ ; Hp.c is another energy scale factor (in addition to (/3c)" 1 ); 

ii) Xp c = (obtained without use of (/3 C ) _1 and a c ), where the preferred direction, 

= ( ^Fr ) ' characterizes thermodynamic continuity in this diagram ; 

L v ' v p J cp 

iii) Sp c = (also obtained without use of (/3 C ) 1 and a c ), where the critical entropy per 
particle Sp, c — 



{ot) 



> is a preferred direction connected to r ) c (measured on the 

CP 

p;T diagram, Fig. 2a), and S c (unmeasured on the Hp]T diagram, Fig. 2b), throughout the 
relation sp. c = -f' c vp. c - S' c , leading to s^ c + x^ c = Y C Z C - j± c , with -j± c = Z c . 



The product Y C Z C = kg cj CP is independent of the reduction process, then is a 

particle property, characteristic of each one-component fluid. De facto, the non-dimensional 
pair {Z C ,Y C Z C } (obtained without use of (/? c ) and a c ) is associated to the two preferred 
directions of $j PiV=cte {jp,v=cte,T,Vp) = 0, where jp,v=cte (T, Vp) = (jf) v=cte is the Grand 
potential per particle of a fluid maintained at constant volume. 

Finally, the extended set {T c ;p c ; Z c ; Y c ; /i^ c ; c\ CP is the complete set from which we 
are able to calculate all the other characteristic critical parameters using linearized thermo- 
dynamics. For example, from U = TS + G + J, we obtain = s Ct p + + ^p^ ^ and 
then u* 3 = n* c - + s* c3 - Z c . 



6 



EUROPHYSICS LETTERS 





P = P C 










(G) 














/ 


— s p. 




/ 






/' 


-l' v \ 
> c p.c n—n r 




temperature ^ j 

Figure 2 - Preferred directions crossing critical point in fip',p (a), \i^\T (b), and p;T (c) diagrams ob- 
tained from projections of the 3D characteristic surface p(T, /ip)and corresponding isothermal (blue), 
isobaric (black), isochoric (green), and isochemical potential (yellow) isoclines. Non- homogenous 
two-phase domain corresponds to the liquid- vapor equilibrium (LVE) (red) line. See also [18]. 



We can now reformulate critical thermodynamics, at the scale of the microscopic interac- 
tion cell of volume v c j — fc -^ Tc filled with -g- particles, where we expect that all the system 
information is contained. Multiplying then the above particle relations by -g-, we obtain 

A* i 

Jc,i — 1 

<j-l = Y c -X c (3) 
<j = s *cj + - 1 

where X c = As expected, Z c disappears in the above equations. X c , as well as Y c , are 

two characteristic properties per CIC. 

Since the Grand potential (which favorizes the local ordering) takes master value at crit- 
icality, j* j = —1, whatever the selected CIC volume, it seems natural to postulate that the 
entropy at criticality (which favorizes the local disorder) can also take a master constant value, 
s*j — const. From Eqs. 0, fx* j — u* j is then also constant. More generally, all the remain- 
ing critical free energies per CIC volume, u* 7 , g*j = H* c j, h* j = u* j — j* j = 1 + u* c 7 , and 
a* j — j* j + g*j = — 1 + fi* j, are also master constants, except a possible common constant 
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value corresponding to an energy translation. Finally, the two dimensionless preferred direc- 
tions X c and Y c are related. The remaining set of the characteristic parameters which contains 
the complete information at the CIC scale is well Q" lm . Our above suggestion complements 
the third law of thermodynamics (s-p — at T — 0), and defines particle entropy and particle 
free energy from (unknown) constants of proportionality to Z c at exact CP, in conformity 
with basic thermodynamic principes. Therefore (/3 C ) _1 is the unique fluid-dependent scale 
factor for energy, as (a c ) _1 is the unique fluid-dependent scale factor for length, as initially 
postulated in [4]. 

As a conclusion, from a formulation of critical thermodynamics in units of the properties of 
the CIC volume, the above analysis shows that the appropriate scale dilatation of the physical 
fields seems adequate to observe scaling of their asymptotic critical singularities (with only one 
adjustable parameter to account for quantum effects). To describe the fluid singular behavior 
within the preasymptotic domain, scale dilatations of two independent fields can then be used 
as controlled simplifications of linear combinations of three fields [19] associated to the density 
formulation of non-symmetrized thermodynamic potentials at finite distance to the CP. 
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